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In order to describe quantum heat engines, here we systematically study isothermal and isochoric processes 
for quantum thermodynamic cycles. Based on these results the quantum versions of both the Carnot heat engine 
and the Otto heat engine are defined without ambiguities. We also study the properties of quantum Carnot and 
Otto heat engines in comparison with their classical counterparts. Relations and mappings between these two 
quantum heat engines are also investigated by considering their respective quantum thermodynamic processes. 
In addition, we discuss the role of Maxwell's demon in quantum thermodynamic cycles. We find that there is 
no violation of the second law, even in the existence of such a demon, when the demon is included correctly as 
part of the working substance of the heat engine. 



PACS numbers: 05.90.+m, 05.70.-a, 03.65.-w, 51.30.+i 

I. INTRODUCTION 

Quantum heat engines (QHEs) 01 [H produce work us- 
ing quantum matter as their working substance. Because of 
the quantum nature of the working substance, QHE have un- 
usual and exotic properties. For example, under some condi- 
tions, QHE can surpass the maximum limit on the amount of 
work done by a classical thermodynamic cycle |13l|4j] and also 
surpass the efficiency of a classical Carnot engine cycle fU. 
QHEs offer good model systems to study the relation between 
thermodynamics and quantum mechanics. Meanwhile, they 
can highlight the difference between classical and quantum 
thermodynamic systems, and help us understand the quantum- 
classical transition problem of thermodynamic processes JfJ. 

The classical Carnot heat engine is a well-known machine 
that produces work through thermodynamic cycles. The ther- 
modynamic properties of the four strokes of each cycle are 
simple and demonstrate the universal physical mechanism of 
heat engines. Current studies |HHSBS!H,[M[H 
[l2l [T3I [Till on QHE mostly focus on the quantum analogue 
of classical Carnot engines, i.e., the quantum Carnot engine 
(QCE). The quantum Otto engine (QOE) is another interest- 
ing case of aQHE, which is also attracting considerable atten- 
tions 0! [3 mini [3. However, we find that there is no 
universal and consistent definition of the QCE and the QOE 
in literatures (see, e.g., Refs. lHl9l[Iol[l9U20ll), and thus the 
properties of QCEs and QOE are not always addressed ade- 
quately and clearly. 

Any QHE cycle consists of several basic quantum ther- 
modynamic processes, such as quantum adiabatic processes 
(which have been clarified in many references, e.g., ll2llD . 
quantum isothermal processes and quantum isochoric pro- 
cesses. This paper begins by clarifying the concepts of 
isothermal processes, isochoric processes and effective tem- 
peratures in their quantum mechanical pictures. Then we sys- 
tematically study the general properties of a quantum ana- 
logue of a Carnot engine. The difference between a QCE and 
its classical counterpart is indicated clearly. We also study the 



QOE based on its basic quantum thermodynamic process and 
analyze the relation between these two types of QHEs. Here 
we assume that the processes of our thermodynamic cycles 
are infinitely slow, i.e., the time interval of each process is as- 
sumed to be very long. Accordingly, the output power is very 
small. This is also the requirement of quasi-static processes. 
Assuming fast cycles would increase the output powers, but 
at the expense of reduced engine efficiency. Also, some ex- 
perimentally realizable physical systems, which can be used 
to implement our QCE and QOE, are discussed. Furthermore, 
based on our generalized QOE model, we demonstrate that 
there is no violation of the second law, even in the presence of 
a Maxwell's demon. 

Our paper is organized as follows: In Sec. II we give 
a clear definition of quantum isothermal and isochoric pro- 
cesses based on the quantum identification of work performed 
and heat exchange. In Sec. Ill we discuss the QCE cycle and 
calculate the work done during this cycle and its operation ef- 
ficiency. In Sec. IV we discuss the QOE cycle and compare 
it with the classical Otto engine cycle. In Sec. V we com- 
pare these two kinds of QHEs and study the relation between 
them. In Sec. VI we give some examples of these two kinds 
of QHEs considering experimentally-realizable physical sys- 
tems. In Sec. VII we discuss the QOE and Maxwell's demon. 
Conclusions and remarks are given in Sec. VIII. 



II. BASIC QUANTUM THERMODYNAMIC PROCESS 
A. Quantum first law of thermodynamics 

To define quantum isothermal and quantum isochoric pro- 
cesses, we need to first consider the working substance. An ar- 
bitrary quantum system with a finite number of energy levels 
is used here as the working substance (see Fig. 1). (Of course, 
this can be generalized to systems with an infinite number of 
energy levels). The Hamiltonian of the working substance can 
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be written as 



H = Y d E n \n) (n\ 



(1) 



where \n) is the n-th eigen state of the system and E n is its 
corresponding eigen energy. Without loss of generality, we 
choose the eigen energy of the ground state |0) as a refer- 
ence point (see Appendix A). Then the Hamiltonian ([1]) can 
be rewritten as 



H = ^2(E n - E ) \n) (n\ 



(2) 



Below we will show that it is convenient for our discussion 
about QHEs to use the Hamiltonian |0. The internal energy 
U of the working substance can be expressed as 



U= (H) = Y,P n E n , 



El 



(3) 



for a given occupation distribution with probabilities P n in the 
?i-th eigen state. 

To clearly define quantum isothermal and isochoric pro- 
cesses, we need to identify the quantum analogues of the heat 
exchange dQ and the work performed &W . From Eq. (O we 
have 



dU = J2 \ E n dPn + Pn dE n ] 



(4) 



In classical thermodynamics, the first law of thermodynamics 
is expressed as 



dU =dQ+dW 1 



(5) 



where dQ = TdS, and dW = J2i Y ^ d Vi t is the tem- 
perature and S is the entropy; yi is the generalized coordinates 
and Yi is the generalized force conjugated to y t . Due to the re- 
lationship S = —fog J2i Pi m -Pi between the entropy 5* and 
the probabilities Pi, we can make the following identification 
MS 



dQ = J2 E ndPn, 

n 

dW = Y, P ndE„. 



(6) 

(7) 



Equation implies that the work performed corresponds to 
the change in the eigen energies E n , and this is in accor- 
dance with the fact that work can only be performed through 
a change in the generalized coordinates of the system, which 
in turn gives rise to a change in the eigen energies [4, 23]. 
Thus the quantum version of the first law of thermodynamics 
dU = dQ + dW just follows from Eq. © with the quan- 
tum identifications of heat exchange and work performed in 
Eqs. © and (0. Different from dQ = TdS, which is appli- 
cable only to the thermal equilibrium case, below we will see 
that Eqs. (|6]l and (Q are applicable to both the thermal equi- 
librium case (see e.g., Eq. d36l l) and the nonequilibrium case 
(see e.g., Eq. (flTTi). 



FIG. 1: (Color online) Schematic diagram of multi-level quantum 
system as the working substance for a QHE. E% and E l n are the n-th 
eigen energy of the working substance in the two isochoric processes. 



B. Quantum isothermal process 

Let us now consider the quantum versions of some thermo- 
dynamic processes. First we study quantum isothermal pro- 
cesses. In quantum isothermal processes, the working sub- 
stance, such as a particle confined in a potential energy well, is 
kept in contact with a heat bath at a constant temperature. The 
particle can perform positive work to the outside, and mean- 
while absorb heat from the bath. Both the energy gaps and 
the occupation probabilities need to change simultaneously, 
so that the system remains in an equilibrium state with the 
heat bath at every instant. Specifically, let us consider a two- 
level system with the excited state |e), the ground state \g), 
and a single energy spacing A. In the quasi-static quantum 
isothermal process, the ratio r = P e /P g of the two occupation 
probabilities, P e and P g , must satisfy the Boltzmann distribu- 
tion r = P e /Pg = exp[— (3A(t)] and also the normalization 
condition P e + P g = 1. A(<) changes slowly with time t, and 
accordingly r can be written as 



"(*) = 



Pn 



(8) 



where j3 — 1/ksT, ks is the Boltzmann constant and T is the 
temperature. In a sufficiently slow process, at every instant the 
system remains in thermodynamic equilibrium with the heat 
bath. 



C. Effective temperature 

We can also define an effective temperature T c s for any 
two-level system according to the ratio r(t) and the level spac- 
ing A(t). For a two-level system with energy levels E e and 
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TABLE I: Quantum versus classical thermodynamic processes. Here we use "INV" to indicate the invariance of a thermodynamic quantity 
and "VAR" to indicate that it varies or changes. U is the internal energy of the working substance; T, P, E n , P n are defined in Sec. II. The 
working substance of the classical thermodynamic processes considered here is the ideal classical gas. 





isothermal process 


isochoric process 


adiabatic process 


Classical 


Heat absorbed or released. 

Work done. 
INV: U,T VAR: P, V 


Heat absorbed or released. 

No work done. 
INV: V VAR: P, T 


No heat exchange. 
Work done. 
VAR: P,T, V 


Quantum 


Heat absorbed or released. 
Work done. 
INV: T VAR: U,E n ,P n 


Heat absorbed or released. 
No work done. 
INV: E n VAR: P n ,T cB 


No heat exchange. 
Work done. 
INV: P n VAR: E n , T e g 



Eg, even in a non equilibrium state, we can imagine that it is 
in a virtual equilibrium state with the effective temperature 



T, 



1 



elf 



A(f) 

fcs/3 eff k B 



■4 



(9) 



as long as the level spacing A(t) and the energy level dis- 
tributions P g and P e are known. Of course, Eq. (O cannot be 
directly generalized to the case with more than two levels. For 
example, for a three-level system with occupation probabili- 
ties P a , Pfc, and P c in three states denoted by \a), \b), and |c), 
if the two level spacings A a b(t) and Ab c {t) do not satisfy the 
relation 



1 



1 



ln Pa = 

A ab (t) n p fc A bc (t) 



Pr 



(10) 



we cannot define a unique effective temperature. The subset 
{ | a) , b) } can have an effective temperature defined by Eq. (O, 
while the subset {\b), c}} would have a different effective 
temperature. We will discuss this point in detail in a QCE 
cycle in Sec. III. 



D. Quantum isochoric process 

A quantum isochoric process has similar properties to that 
of a classical isochoric processes. In a quantum isochoric pro- 
cess, the working substance is placed in contact with a heat 
bath. No work is done in this process while heat is exchanged 
between the working substance and the heat bath. This is the 
same as that in a classical isothermal process. In a quantum 
isochoric process the occupation probabilities P n and thus the 
entropy S change, until the working substance finally reaches 
thermal equilibrium with the heat bath. In classical isochoric 
process the pressure P and the temperature T change, and the 
working substance reaches thermal equilibrium with the heat 
bath only at the end of this process. For example, if the work- 
ing substance is chosen to be a particle confined in a infinite 
square well potential, no work is done during a quantum iso- 
choric process when heat is absorbed or released, and the oc- 
cupation probabilities in every eigen state satisfy Boltzmann 
distribution at the end of the isochoric process. 



E. Quantum adiabatic process 

A classical adiabatic thermodynamic process can be formu- 
lated in terms of a microscopic quantum adiabatic thermody- 
namic process. Because quantum adiabatic processes proceed 
slow enough such that the generic quantum adiabatic condi- 
tion is satisfied, then the population distributions remain un- 
changed, dP n = 0. According to Eq. (|6), dQ = 0, there is 
no heat exchange in a quantum adiabatic process, but work 
can still be nonzero according to Eq. (0. A classical adiabatic 
process, however, does not necessarily require the occupation 
probabilities to be kept invariant. For example, when the pro- 
cess proceeds very fast, and the quantum adiabatic condition 
is not satisfied, internal excitations will likely occur, but there 
is no heat exchange between the working substance and the 
external heat bath. This thermodynamic process is classical 
adiabatic but not quantum adiabatic. Thus it can be verified 
that a classical adiabatic process includes, as a subset, a quan- 
tum adiabatic process; but the inverse is not valid lfl7tl . 

The properties of both classical and quantum thermody- 
namic process are listed in Table I, to facilitate the comparison 
between these two processes. The table indicates if heat is ab- 
sorbed or released (first row), if work is done (second row), 
and which quantity varies (indicated by "VAR") and which 
are invariant (indicated by "INV") in the third row. 



III. QUANTUM CARNOT ENGINE CYCLE 

In the previous section, we defined quantum isothermal pro- 
cesses. Based on this definition, in this section, we study the 
QCE cycle and its properties. The QCE cycle (see Fig. 2 for an 
example of a QCE based on a two-level system), just like its 
classical counterpart, consists of two quantum isothermal pro- 
cesses (A — > B and C — > D) and two quantum adiabatic 
processes (B — > C and D — > A). During the isothermal 
expansion process from A to B, the particle confined in the 
potential well is kept in contact with a heat bath at tempera- 
ture Th, while the energy levels of the system change much 
slower than the relaxation of the system, so that the particle is 
always kept in thermal equilibrium with the heat bath. Below, 
we consider both cases: two-level and multi-level systems. 
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(A) quantum Carnot engine 
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(C) T-S diagram of QCE (CCE) 
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FIG. 2: (A): A schematic diagram of a quantum Carnot engine based 
on a two-level quantum system. A is the level spacing between the 
two energy levels. P e is the occupation probability in the excited 
state. The process from A to B (C to D) is the isothermal expan- 
sion (compression) process, in which the working substance is put 
in contact with the high (low) temperature heat bath. The processes 
from B to C and from D to A are two adiabatic processes. (B): 
Pressure-Volume (PV) diagram for a classical Carnot engine with 
ideal gas as the working substance. The process from 1 to 2 (from 
3 to 4) is the classical isothermal expansion (compression) process 
with temperature TJ, (T;), and the process from 2 to 3 (4 to 1) is the 
classical adiabatic expansion (compression) process. V% and Vz are 
the volume of the working substance at 2 and 3 respectively. (C): 
Temperature-Entropy (T — S) diagram [27] for both for a quantum 
Carnot engine based on a two-level quantum system and a classical 
Carnot engine with ideal gas as the working substance. This T — S 
diagram bridges the quantum and classical Carnot engine. 



A. Thermodynamic reversibility of the quantum Carnot 
engine cycle 

It is well known that quantum mechanical reversibility are 
associated with quantum mechanical unitary evolution. Dif- 
ferent from quantum mechanical reversibility, thermodynamic 
reversibility accompanies the heat bath and the effective tem- 
perature of the working substance. In this paper we focus on 
the thermodynamic reversibility. 

We emphasize that, in order to ensure that the cycle is 
thermodynamically reversible, two conditions on the quantum 
adiabatic process are required: (1) after the quantum adiabatic 
process (B — > C), we can use an effective temperature 7} 
to characterize the working substance, i.e., the working sub- 
stance still satisfies the Boltzmann distribution after the quan- 
tum adiabatic process; and (2) the effective temperature Ti of 
the working substance, after the quantum adiabatic process, 
equals the temperature Ti of the heat bath of the following 
quantum isothermal process (C — ► D). When either condi- 
tion is not satisfied, a thermalization process $2a\ l25ll of the 
working substance is inevitable before the quantum isother- 
mal process (C — ► D). In the thermalization process, the 
total entropy increase of the working substance plus the bath 
is nonzero. Hence, this thermalization process is irreversible. 

It can be proved that the above two conditions are equiv- 
alent to the following two conditions: (i) all energy gaps 
were changed by the same ratio in the quantum adiabatic 
process, i.e., E n (B) - E m {B) = \[E n {C) - E m {C)], and 
E n (A) - E m (A) = X[E n (D) - E n (D)}, (n = 0, 1, 2, • • • ); 
and (ii) the ratio of the change of the energy gaps in the adia- 
batic process (from BtoC or from A to D) must equal the ra- 
tio of the two temperatures of the heat baths, i.e., A = Th/Ti. 

First, it is easy to see that these two conditions, (i) and (ii), 
listed right above are sufficient for the previous conditions ( 1 ) 
and (2) presented initially. Next, we prove that the two con- 
ditions (i) and (ii) are also necessary for the two conditions 
(1) and (2). Let us assume that the working substance is in 
equilibrium with a heat bath at temperature 7), at the instant 
B before the adiabatic process (B — > C). In this case the 
quantum state is described by a density operator 

P( B ) = \ E eX PH^n(5)] \n{B)) (n(B)\ . (11) 

n 

After the adiabatic process is completed, at instant C in Fig. 2, 
the eigen energies of the working substance become E n {C), 
and the working substance reaches an effective temperature 
Ti. The occupation probabilities B n of the working substance 
are kept unchanged during the adiabatic process (B — > C), 
and they satisfy the Boltzmann distribution. Thus, for any 
eigen states \n) and m), the occupation probabilities P n and 
P m satisfy 

P n (B) = exp[-!3 h E n (B)] 
P m {B) exp[-(3 h E m (B)} 

_ P n (C) = cxp[-/^ n (C)] 

P m (C) exp[-/3,£ m (C)r 1 ; 
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That is, 



E n {C) - E m (C) = ^-[E n (B) - E m (B)], 

J-h 



(13) 



for any m, n. Eq. ( fT3l is just a combination of conditions (i) 
and (ii). Thus we have proved that the two conditions (i) and 
(ii) are sufficient for the previous conditions (1) and (2). 

Hence, we have proven that (i) all energy gaps change by 
the same ratio in quantum adiabatic process, and (ii) this ratio 
equals to the ratio of the temperatures of the two heat baths, 
summarized in Eq. (1131 1. are sufficient and necessary condi- 
tions for the QCE to be thermodynamically reversible. We 
would like to mention that these two conditions (i) and (ii) 
(mathematical requirements) can be satisfied in some realistic 
physical systems. Examples of QCEs based on some concrete 
physical systems will be discussed in Sec. VI. 



B. Work and efficiency of a quantum Carnot engine cycle 

Now we analyze the operation efficiency rj c of the QCE in- 
troduced above. For simplicity, instead of applying Eq. (O, 
we use dQ = TdS to calculate the heat exchange dQ in any 
Quantum Iso-Thermal (QIT) process. Because the tempera- 
ture of the heat bath is kept invariant in the quantum isother- 
mal process, the heat absorbed Q^ IT and released Q®ut~ m 
the quantum isothermal expansion and compression processes 
can be calculated as follows 



g IT = T h [S(B)-S(A)}>0, 
T l [S(C)-S(D)}>0, 



} QIT = 
i out 



(14) 
(15) 



where Th and T; are the temperatures of the two different heat 
baths, and 



S(i) 



\ - exp[-&g„(z)] 



Z(i) 



(16) 

are the entropies of the working substance at different in- 
stants i = A, B, C, D (see Fig. 3(A)). Here, (3 A B = 
X/ksTh, fic d — 1/ksTi. In obtaining the above result, 
we have used the Boltzmann distribution of thermal equilib- 
rium state, i.e., p = (l/Z)J2 n ex P(—(3E n ) \n) (n|, where 
Z =Trexp(— (3H) is the partition function. Of course, Q^ IT 

and Q^^t can a l so t> e obtained through Eq. © in a quan- 
tum manner (this will become clear later on, see Eqs. ( f36l > and 
(|37||). These equivalent approaches can describe the micro- 
scopic mechanism of a classical Carnot engine cycle. 

Now, we would like to calculate the work Wc done during 
a QCE cycle and its operation efficiency t] c . From Eqs. (fl4] i 
and (IT~5b and the first law of thermodynamics we obtain the 
net work done during a QCE cycle 



W C = Q? n lL - = (T h - Ti)[S(B) - S(A)l (17) 

where we have used the relations S(B) = S(C) and S(A) = 
S(D). This equivalence is due to the fact that the occupa- 
tion probabilities and thus the entropy remain invariant in any 



quantum adiabatic process. The efficiency r) c of the QCE is 



'1c 



Q 



(18) 



which is just the efficiency of a classical Carnot engine. From 
Eq. ( fT3l we see that the ratio of the temperature in the effi- 
ciency (Eq. (fT8l ) of the QCE can also be replaced by the ratio 
of the energy gaps 



Vc 



E n {C) — E m (C) 
E n (B) — E m (B) ' 



(19) 



This expression of the efficiency rj c in terms of the ratio of the 
energy ga ps r esembles that of a QOE for a multi-level case in 
Refs. yl ll7[] (see also Eq. d24i > below). However, in spite 
of the apparent similarities between these two expressions for 
the efficiencies, we emphasize that they are quite different. 
Here, E n (B) - E m (B) and E n (C) - E m (C), in Eq. tO, 
are the energy gaps at the beginning (B) and at the end (C) 
of the quantum adiabatic expansion process (B — ► C). In 
the expression for the efficiency r] for a multi-level QOE, 
however, the energy gaps are those in two quantum isochoric 
processes 1 3l 1 1711 . Hence, the efficiency in Eq. ( fT9l ) for a QCE 
is quite different from that for a QOE, even though they both 
look similar. Further discussions on this will be given in Sec. 
IV. 

In order to extract positive work from the bath, Eq. ( fT8l im- 
poses a constraint, Th > Ti, on the temperatures of the two 
heat baths. This constraint, known as the positive-work con- 
dition (PWC), is the same as that of its classical counterpart. 
What is more, the schematic temperature-entropy (T — S) di- 
agrams for both a QCE cycle and a classical Carnot engine 
cycle are the same (see Fig. 2(C)). For the above reasons, we 
believe it is convincible that our QCE model is a quantum me- 
chanical analogue of a classical Carnot engine. We compare 
the properties of a QCE and a classical Carnot engine and list 
them in Table II. 



C. Internal energy 

It is well known that an ideal classical Carnot engine cycle 
consists of two classical isothermal and two classical adiabatic 
processes. When the working substance is the ideal gas, the 
internal energy of the working substance remains invariant in 
the classical isothermal process, because the internal energy 
of the ideal gas depends on the temperature only. This as- 
sumption for classical isothermal processes based on classical 
ideal gas could be true for a classical Carnot engine using a 
working substance other than an ideal gas. But in the quan- 
tum version, the quantum isothermal and quantum adiabatic 
processes should be redefined microscopically based on quan- 
tum mechanics. In principle, the classical result could change 
when considering the quantum nature (discrete energy levels) 
of the working substance. 

We now would like to verify whether the internal energy of 
the working substance remains invariant during the isothermal 
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TABLE II: Quantum Carnot engines versus classical Carnot engines. Here "CIT" refers to "Classical Iso-Thermal process" while "CA" is an 
abbreviation for a "Classical Adiabatic process". "QIT" and "QA" refer to "Quantum Iso-Thermal process" and "Quantum Adiabatic process", 
respectively. V2, V3, E n (B), and E n (C) are defined in Fig. 2; 7 is the adiabatic exponent 12711 . 





strokes 


requirement on the CA and QA 


efficiency 


positive-work condition 


Classical 


CIT-CA-CIT-CA 


T, _ rV 2 i 7 -i 


u = i-£ 


T h > T ; 


Quantum 


QIT-QA-QIT-QA 


TV E n (C)-E m (C) r u 

T h ~ E n (B)~E m (B) forVm,n 




T h > Ti 



process. At the four instants A, B, C and D of the QCE cycle 
(see Fig. 2), the internal energies are respectively 



U(i) = Tr[p(i)H(i)l i = A,B,C,D. 



(20) 



In appendix B we prove that U(A) ^ U(B) and U(C) ^ 
U{D) for some QCE cycles based on several experimentally 
realizable systems. Hence, in the quantum version of a Carnot 
engine, we cannot simply assume that the heat absorbed (re- 
leased) by the working substance equals to the work done by 
(on) the working substance in the isothermal process, as we 
do in classical Carnot engine with the ideal gas as the work- 
ing substance. This observation is crucial for the following 
discussion. 

Here we would like to indicate that, the quantum isoener- 
getic process in Refs. (JtL [^1 is not a quantum analogue of the 
classical isothermal process of classical Carnot engine, be- 
cause it requires the temperature of the heat bath to change. 
Thus the thermodynamic cycle described in Refs. [Ml is 
actually not a QCE cycle. 



IV. QUANTUM OTTO ENGINE CYCLE 

In practice, the heat engines most widely used in auto- 
mobiles, the internal combustion engine, operate using Otto- 
cycle engines IU8I1 . which consist of two classical isochoric 
and two classical adiabatic processes. Similar to the Carnot 
engine, the quantum analogue of the classical Otto engine is 
also proposed in Refs. (SH The QOE cycle 

consists of t wo q uantum isochoric and two quantum adiabatic 
processes ||4j, |l5l |l6|, ll7[] (see Fig. 3 for a schematic diagram 
of QOE based on a two-level system). 



A. Work and efficiency 

In the quantum isochoric heating process from A to B (see 
Fig. 3), no work is done, but heat is absorbed. The heat Q^ IC 
absorbed by the working substance is 

QT = E C E " dP » = E E n [Pn(B)-P n (A)], (21) 

n jA 

where is the n-th eigen energy of the system in the quan- 
tum isochoric heating process from A to B. Similarly, we 



obtain the heat released to the low temperature entropy sink 
in the quantum isochoric cooling process from C to D 



E n dP n = E K [Pn(C) - Pn(D)}, 



(22) 

where E l n is the n-th eigen energy of the system in the quan- 
tum isochoric cooling process. We would like to point out that 
in calculating Q^ IC and Q^ , we cannot apply Eqs. ( TBI ) 
and ( [15] ), because &Q = TdS is only applicable to thermal 
equilibrium case, while in the quantum isochoric process, the 
heat bath and the working substance are not always in ther- 
mal equilibrium, i.e., this process is not thermodynamically 
reversible (for a detailed discussion see below). 

As mentioned above, in order to construct a QCE, all energy 
gaps must be changed by the same ratio in quantum adiabatic 
process. But for a multi-level QOE, there is no such a con- 
straint (see Ref. [17]) because we do not have to ensure the 
reversibility of the QOE cycle. Nevertheless, to compare the 
QOE with the QCE, we only consider a special case of QOE 
where all its energy gaps change by the same ratios as in the 
quantum adiabatic processes, i.e., E„ — E' r l n — a(E l n — E l m ), 
(n = 0,1,2,.. .). When we choose £# = E l = 0, i.e., the 
ground state eigen energies as the energy reference point, we 
have Ell = a ^n- Similar to the QCE, the occupation distribu- 
tion remains invariant in the two quantum adiabatic processes, 
i.e., P n (B) = P„(C) and P n (A) = P n (D), and accordingly 
the entropy remains invariant in the quantum adiabatic pro- 
cesses S(B) = S(C) and S(A) = S(D). 

Based on this fact, and Eqs. d2"TT i and i22i . we obtain the net 
work Wq done during a QOE cycle 



Wo 



)Qic 

i out 



(E*-E l n )[P n (B)-P n (A)}. 

(23) 



and the operation efficiency ?/ of the QOE cycle 



'to 



Wo 

Q 



QIC 



= 1 - 



E h - E h 



= 1 - 



1 



(24) 



Here, a > 1 since E% > E l n . This result, which stands for a 
special multi-level QOE (all energy gaps change by the same 
ratio in the quantum adiabatic process), is a generalization of 
the two-level QOE dflj. Let us recall the PWC of the spe- 
cial multi-level QOE [17] mentioned above. From Eq. d23l ). 
the PWC for the special multi-level QOE is 



T h > a Tj. 



(25) 
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(A) quantum Otto engine 

A 

V-X(B)=T h 
-T(C) 




Pi 

e 



Ph p 

e e 



(B) classical Otto engine 




h 1 
(C) T-S diagram of QOE (COE) 



B (2) 



A (1) 




D 1(4) 



S (A ) S (B ) 



FIG. 3: (Color online) (A): A schematic diagram of a quantum Otto 
engine based on a two-level quantum system. Dotted red and dashed 
black curves refer to the isothermal processes. and A; are the 
level spacings of the two-level system during the two isochoric pro- 
cesses (A — > B and C — > D). P e and Pj are the occupation 
probabilities in the excited state. The processes from A to B and 
from C to D are isochoric processes, while the processes from B to 
C and from D to A are quantum adiabatic process. T(i), (i = A, 
B, C, D) are the effective temperatures of the working substance 
at instant i, and T(B) = T h , T(D) = T h (B): Pressure- Volume, 
PV, diagram of a classical Otto engine. Vj» and Vi are the volume of 
the working substance in the two classical isochoric processes. The 
process from 1 to 2 (3 to 4) is a classical isochoric heating (cooling) 
process, and the process from 2 to 3 (4 to 1) is the classical adiabatic 
expansion (compression) process. The temperature of the working 
substance at 2 and 4 equal to the temperatures of the two heat baths, 
Th and Tj, respectively. (C): Temperature-Entropy (T — S) diagram 
12711 for both for a quantum Otto engine based on a two-level quan- 
tum system and a classical Otto engine with ideal gas as the working 
substance. This T — S diagram bridges the quantum and classical 
Otto engine. 



This is obviously different from that of a QCE, where the 
PWC is simply T h >T\. 

Note that the first QHE model, initially proposed in Ref. 
Ql], is actually a QOE, because its efficiency and its PWC are 
given by r\ = 1 — v\jv v and T\ > {v v Iv\)Tq, where v\ and 
v v are the two energy gaps of the working substance, and T\ 
and Tq are the temperatures of the two heat baths, respectively. 



B. Classical versus quantum Otto engines 

Below we prove that the operation efficiency rj Q in Eq. (f24b 
of a QOE, also equals to the efficiency ?7q L of a classical Otto 
engine. For simplicity, here we only consider a two-level sys- 
tem as the working substance (the result can be generalized to 
multi-level systems if all the eigen energies of the multi-level 
system change in the same ratios as in the quantum adiabatic 
process lr26Tl ). For a two-level system (see Fig. 3), when the 
temperature T of the heat bath is fixed, the occupation prob- 
ability P e of the excited state |e) in thermal equilibrium is a 
monotonically decreasing function of the level spacing A be- 
tween the two energy levels 0, EH H] ■ Its inverse function 
reads 



Ag = k B Tg In 



PS 



h,l. 



(26) 



As mentioned above, the efficiency rj Q d24i l of a QOE cycle 
represented by the rectangle (A-B-C-D) (see Fig. 3) is 



77o 



A h 



(27) 



From Eqs. d26| i and (IZTt and Fig. 3, we can see that the effi- 
ciency rj of the QOE cycle can be rewritten as 



Vo = 1 - 



T(C) 
T(B) 



= 1 - 



T(D) 
T(A) ' 



(28) 



where T(i), with i — A, B, C and D, are the effective tem- 
peratures of the working substance at the instants A, B, C 
and D indicated in Fig. 3. Here, we have used the relation 
A(C)/A(B) = T(C)/T(B) because of the fact P e (C) = 
P e (B) = P£ and Eq. ([35}, and similarly A(D)/A(A) = 
T(D)/T{A). In the QOE cycle, the effective temperatures 
T{B) and T(D) of the working substance at instants B and 
D equal the temperatures of the two heat baths T(B) = Th, 
T(D) = T,. 

As for a classical Otto engine, the classical efficiency t]q L 
is 1271 



(29) 



where Vi and Vh are the volumes of the classical ideal gas 
in the two classical isochoric processes (see Fig. 3), and 7 
is the classical adiabatic exponent [27]. Because TV 1 ' 1 is 
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TABLE III: Quantum Otto engine versus classical Otto engine. Here "CIC" and "CA" refer to "Classical Iso-Choric processes" and "Classical 
Adiabatic processes", respectively; "QIC" and "QA" refer to "Quantum Iso-Choric processes" and "Quantum Adiabatic processes", respec- 
tively. Also, Vh and Vi are the volumes of the working substance (classical ideal gas) in the two classical isochoric processes; 7 is the classical 
adiabatic exponent ['27]. T(i) (i = A, B, C, D) and T k (k = 1, 2, 3, 4) are defined in Fig. 3 





strokes 


efficiency 


positive-work condition 


Classical 


CIC-CA-CIC-CA 


CL _ (V lL ^ 1 -l_ , T(C) _ 1 T(D) 
'lO — 1 I V, I — 1 T[B) — 1 T{A) 




Quantum 


QIC-QA-QIC-QA 


n — 1 — Ai — 1 Zi — 1 T -4 


T h >T l {%) 



constant during a classical adiabatic process, we can therefore 
eliminate the volumes in Eq. (|29i l as 



T7 CL 



T 2 



(30) 



where T\, T 2 , T 3 and T 4 are the temperatures of the working 
substance at instants 1, 2, 3, and 4, and the temperatures at 
instants 2 and 4 equal to the temperatures of the two heat baths 
T 2 = T 4 = T\. This is in very good agreement with 
the result of a QOE cycle in Eq. (f28t . Thus we proved that 
the efficiency of a QOE also equals its classical counterpart. 
Moreover, similar to Carnot engines, we plot the schematic 
temperature-entropy (T — S) diagrams for both a QOE cycle 
and a classical Otto engine cycle in Fig. 3(C). The similarity 
of the two T — S diagrams also support our definition of the 
QOE. Comparisons between the QOE and the classical Otto 
engine are listed in Table III. 



C. An alternative quantum Otto engine 

Before concluding this section, we would like to revisit an 
alternative QOE cycle similar to that in Ref. lf28tl and that 
given most recently by us 1 18]. As illustrated in Fig. 4, we 
consider two two-level systems (two qubits) as the working 
substance, which are denoted by qubit S and qubit D, and 
the level spacings of the two qubits are Ag and A^, respec- 
tively. The temperatures of the two heat baths are Tg and Tjj. 
Without loss of generality, here we choose T5 > To and 
As > A D . 

The alternative QOE cycle consists of two steps: 1) let the 
two qubits decouple from each other and contact their own 
heat baths respectively until they reach thermal equilibrium 
with these two heat baths; and 2) switch on the interaction 
between the two qubits and implement a SWAP operation be- 
tween them. These two steps are shown in Fig. 4(a) and 4(b). 
The density matrices of the two qubits after step 1 are 

f\0) = y. [|0) « (0| > + ex p[-^ A *l I 1 ), < 1 L-] . (* = s > D )> 

(31) 

where j3 i — l/hsT, and for simplicity we have chosen the 
eigen energy of the ground state as a reference point. After 




FIG. 4: (Color online) A schematic diagram illustrating an alterna- 
tive quantum Otto engine. This QOE cycle consists of two steps: 
a SWAP between the two two-level systems (two qubits), shown in 
(b), and a thermalization with their respective heat baths, shown in 
(a). Step (c) indicates the transition from (b) to (a). The SWAP oper- 
ation in (b) replaces the two quantum isochoric processes in the QOE 
cycle mentioned in the text and shown in Fig. 3(A). 



step 2, the density matrices become 



Ps(2) 



Z D 



[|0) s (0| s + cxp[- / 3 Z3 A D ]|l) s (l| s ], 



Pd(2) = —[\O) D (O\ D + eM-0s&s] I1M1UP2) 

After these two steps, the two qubits are decoupled and put 
into contact with their own bath, and a new cycle starts. 

The key point of this alternative QOE cycle is that the 
SWAP operation takes place of the two quantum adiabatic 
processes, while the thermolization process takes place of two 
quantum isochoric processes. We can now calculate the heat 
absorbed by the qubit S and heat released by the qubit D in 
step 1: 

Qin = Tr[H Ps (l)]-Tr[Hp s (2)) (33) 
= A s — cxp[-/3 s A s ] - — cxp[-/3 D A D ] 
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Q out - Tr[Hp D (2)]-Ti[Hp D (l 
= A D 



(34) 



— exp[-/3 s A s ] - —exp[-/3 D A D ] 



The operation efficiency r\ of the QHE cycle can be calculated 
straightforwardly 



_ Qin ~ Qout _ A/) 



(35) 



and the PWC is T s > (A S /A D )T D . Thus, this two-step 
cycle is an alternative QOE cycle. We will revisit, in more 
detail, this alternative QOE cycle in Section VII. 

Let us here mention an alternative QCE of two qubits based 
on a SWAP operation. This alternative QCE cycle consists of 
three steps: 1) let the two qubits decouple and contact their 
own heat baths and both experience quantum isothermal pro- 
cesses, 2) switch on the interaction between the two qubits and 
implement a SWAP operation, and 3) let the two qubits de- 
couple from each other and also decouple the two qubits from 
their heat baths, and subject them to a quantum adiabatic pro- 
cess. A similar QOE cycle and a QCE cycle have been studied 
in Ref. [28], where the SWAP operation was decomposed into 
three CNOT operations. 




P! 



p': p. 



FIG. 5: A Quantum Carnot engine cycle can be modeled as infi- 
nite number of small quantum Otto engine. Here, Pg, P l e ,A, B, C, 
D, and T%, T\ are defined in Fig. 2(A). The small rectangles inside 
the QCE cycle (A-B-C-D) represent small QOE cycles. The tem- 
peratures of the two heat baths of these QOE cycles are Th and TJ, 
respectively, which are the same as that of the QCE cycle. Similar 
discussions see also Refs. 



V. RELATIONS BETWEEN QUANTUM OTTO AND 
QUANTUM CARNOT CYCLES 

In this section we discuss the relation between a quantum 
Otto engine cycle and a quantum Carnot engine cycle. For 
simplicity, here we use a two-level system as an example of 
working substance. Our results about QHEs based on a two- 
level system can be generalized to multi-level systems if all 
the eigen energies of the multi-level system change by the 
same ratios in the quantum adiabatic processes B26I1 . 



A. Quantum Carnot cycle derived from quantum Otto cycles 



From Eqs. ( f36b and ( f38l > we can see that the efficiency of an 
infinite number of infinitesimal QOE cycles is 



II" 



V = 



= 1 - 



II 



(39) 



This is the efficiency of the QCE in Eq. ([T8l >. Thus we have 
proved that a QCE cycle can be modeled as an infinite num- 
ber of infinitesimal QOE cycles. However, a finite QCE cycle 
and a finite QOE cycle cannot be equivalent because one is 
reversible and the other one is not. When the two cycles be- 
comes infinitesimal, they can be infinitesimally close to each 
other. 



A QCE cycle can be decomposed into an infinite number of 
small QOE cycles (see Fig. 5) H H H [H. Now we give 
a concise and explicit proof about this observation. The heat 
absorbed and released in the infinite number of infinitesimal 
QOE cycles can be integrated by applying Eqs. (1211 1. ( 122b and 

Qin = Th / lnfP" 1 -!) dP e , (36) 
J PI 

Qout - V / ^(P- 1 -!) dP e . (37) 
J Pi 

Then we obtain the positive work W done during the infinite 
number of infinitesimal QOE cycles {A-B-C-D) by making 
use of Eq. d23l 

W = Qi„ -Qout (38) 

= (T h -Ti) / ln^ 1 -1) dP e . 
J PI 



B. Quantum Otto cycle derived from quantum Carnot cycles 

Conversely, a QOE cycle can also be modelled as an infinite 
number of QCE cycles (see Fig. 6), but the temperatures of the 
two heat baths of these small QCE are different. The quantum 
isochoric process A' — ► B' (C — ► D') can be modelled as 
many small quantum isothermal processes with temperatures 
Tl (T,), Tl (T}\ Tl (TlX . . .(see Fig. 6) 

Tl < Tl < T h 3 < ■ • • < Xf < T h , (40) 
Ti < T; 1 < Tf < ■ • • < T^ 1 < T t N . (41) 

The heat absorbed and released in the infinite number of in- 
finitesimal QCE cycles can be obtained by applying cLQ — 
TdS and Eq. d26j 

J P i c k B In (P e - 1) J P i 
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A 

A; 

A, 



T 1 T 2 



D' 



5 c 




^1 



T 

± 7. 



T 



a 



p! 



p 



p. 



A 
A 



A A 

7X4') 




i 1 i 



c 



B' 



a 



-T(C) 



p! 



p 



FIG. 6: A Quantum Otto engine can be modeled as an infinite num- 
ber of small quantum Carnot engine cycles. Here P e \ P\, T h and Ti 
are defined in Fig. 3(A). The small cycles inside the QOE cycle A 1 - 
B'-C'-D' represent small quantum Carnot cycles. The temperatures 
of the two heat baths (e.g., T h and T; 1 ) of these small QCE cycles 
are different from that (Th and Ti) of the QCE cycle. 



FIG. 7: (Color online) Schematic diagrams for a quantum Carnot 
engine cycle (AB'CD 1 red continuous line) and a quantum Otto en- 
gine cycle (A' B' C' D' blue dashed line) based on a two-level sys- 
tem. Here, points B' and D' are the same for both cycles. The 
temperatures Th and T; of the two heat baths of the two QHE cycles 
and the occupation probabilities and Pi in two quantum adia- 
batic processes are the same. The net work done during a cycle is 
proportional to the area enclosed by the four curves which represent 
the cycles. 



Qo 



k B h\{Pe - 1 



-dS 



dP e , (43) 



where we have used the relation dS = — fcs[lnP e — ln(l — 
P e )]dP e - We then obtain the positive work W done during the 
infinite number of infinitesimal QCE cycles (A'-B'-C'-D' in 
Fig. 6) by making use of Eq. ( TPTI i 



W = Q in - Qout = (A h - A,)(P e ' 1 - P l e ). 



(44) 



From Eqs. d42l > and ( l44b we can see that the efficiency of an 
infinite number of infinitesimal QCE cycle is 



W 

Qin 



= l 



(45) 



This is the efficiency of the QOE in Eq. ( f27b . Thus we have 
proved that a QOE cycle can be modeled as infinite number of 
infinitesimal QCE cycles. 

We would like to mention that the formula dQ = TdS is 
applicable in Eq. d42l > is due to the fact that these infinite num- 
ber of infinitesimal QCE cycles (with heat baths temperatures 
T\, TJ 1 ; T%, T; 2 ; • • • ) are thermal equilibrium (reversible) pro- 
cesses since the entropy increase vanishes during these cycles 
(see Appendix C). However, if we take the QOE cycle as a 
whole, and the heat baths temperatures are Th and 1), this 
process is a nonequilibrium (irreversible) process. Hence, the 
formula dQ = TdS is not applicable here. This is why we 
cannot apply Eqs. ( TBI and (TT3T > in calculating Q^ IC and 
in the QOE cycle in Sec. IVA. 



C. Comparison of work and efficiency for quantum Otto and 
quantum Carnot cycles 

Having clarified the properties of the QCE and QOE, here 
we now compare the thermodynamic properties of a QCE cy- 
cle with that of a QOE cycle and study the relation between 
them. For the two QHE cycles (QCE cycle and QOE cycle), 
we consider the case with the same heat baths (at high and 
low temperatures Th and T/, respectively), and the same oc- 
cupation probabilities (P e /l and P l e , respectively) in the two 
quantum adiabatic processes (see Fig. 7). 

First, let us compare the amounts of positive work done 
in a QCE cycle and a QOE cycle under the same conditions 
defined above. The positive work done during the QCE cycle 
(here denoted by A-B'-C-D' in Fig. 7) can also be written as 



W c 



P! 



[A(T hl P a )-A(T u P e )]dP e 



(46) 



while the work done during a QOE cycle is (here denoted by 
A'-B'-C'-D' in Fig. 7) 



Wo = I (A h - A ; ) dP e 
I Pi 



(47) 



From Fig. 7 we know for any P e £ [P l e ,P^], we have 
A(T h , P e ) - A(T t ,P e ) > A h - A/. Hence, from Eqs. © 
and d4Tb we have Wc > Wo', i.e., under the same conditions, 
the work done during a QCE cycle is more than that during a 
QOE cycle. 

Next we consider the efficiencies of the QCE cycle and the 
QOE cycle under the same conditions. From Eq. ( f27b we 
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know that the efficiency of a small QOE cycle can be rewritten 

as 



'to 



1 - 



Vc 



ksmniiPj)- 1 - 1] 
k B T h hx[(Pl)-i-l] 
d 
' dP e 

IlJL 

T h dP e 



(48) 



1 



111 ln I pi 



(Pe - Pi)- 



It can be verified that the second term on the RHS of Eq. 
is negative. Thus we have proved the inequality 



(49) 



for every small cycle. We conclude that, under the same con- 
ditions, the QCE is more efficient that the QOE, even for any 
finite cycle. 



VI. ILLUSTRATIONS OF QUANTUM CARNOT ENGINE 
AND QUANTUM OTTO ENGINE 

As mentioned above, to construct a multi-level (including 
two-level) Carnot-type QHE, two preconditions (mathemati- 
cal results) are required: (1) all energy gaps change by the 
same ratio in the quantum adiabatic process (when the work- 
ing substance performs work); (2) the ratio of the energy gap 
changes in the quantum adiabatic process should equal the ra- 
tio of the temperatures of the two heat baths, so that the ther- 
modynamic cycle is reversible. Physically, these two condi- 
tions can always be satisfied for a QHE based on a two-level 
system, because there is only one energy gap in a two-level 
system, we can always find a proper effective temperature 
to characterize the working substance. Besides the two-level 
system, the harmonic oscillator and a particle confined in an 
infinite square well potential are two other examples that can 
illustrate the basic properties of the QCE and the QOE. This 
is because in both cases all energy gaps change by the same 
ratio when changing the parameters of the system, and we can 
always use a proper effective temperature to characterize the 
working substance in the quantum adiabatic process, too. Be- 
low we calculate the amount of positive work done during a 
thermodynamic cycle using those working substances. 



A. Two-level systems 

Let us now consider a QHE based on a two-level system, 
e.g., a spin- 1/2 system in an external magnetic field pointing 
along the +z direction, the Hamiltonian of the working sub- 
stance is 



ff T Ls(i) = -MS(i)(|t)(t|-|4)ai), 



(50) 



where B(i) is the strength of the external field at instant i, i = 
A, B, C, D (see Fig. 2 and Fig. 3); M = eh/2mc is the Bohr 
magnon; | J.) and |f) indicate the spin down (excited) state and 



spin up (ground) state respectively. The thermal equilibrium 
state can be written as 



1 



[exp&ATBC*)] |T) <T| 
+ e X p[-/3 i MB(i)]|i>U|], 



(51) 



where Z(i) = exp[/3jMB(i)] + exp[-/3 i MB(i)] is the par- 
tition function at instant i. Applying Eq. ( fT&b we obtain the 
entropy of the working substance 

S TLS (i) = k B H z (i)] + k B P l MB(i) tanh[/3 l M J B(i)]. 

(52) 

Then from Eqs. ( TT4l i, (Q3) and (T% we obtain the net work 
done by the working substance during a QCE cycle 

^c LS - (Th - Ti) [S ThS (B) - S ThS (A)} . (53) 

The net work done during the QOE cycle can be calculated by 
applying Eq. d23l 



W, 



TLS 



(A h - A* 



1 + exp\p h A h ] l + cxp[/3 z A ; 



(54) 



Another example of QHE based on a two-level system is a 
photon-Carnot engine 0, [IH . After performing a similar 
calculation we recover the operation efficiency 1 — 7] /Th in 
Eq. (O and the PWC T h > Tj. Hence, this photon-Carnot 
engine is actually a two-level QCE. 



B. Harmonic oscillator 



For a QHE based on a harmonic oscillator with the eigen 
energies E n (i) = (n — 1/2)Hlj(i), by applying Eq. ([To*l ). we 
derive the entropy S (i) of the working substance as 



5 HU («) = -* fl ln{l-e3rp[-/3 i M*)]} 



(55) 



+k B f3 l Huj(i)- 



1 



exp[/3jMi)] 
and the work done during a QCE cycle as 



1' 



W%° - (T h - T) [S HO (B) ~ S HO (A)] 



(56) 



The net work done during the QOE cycle can be obtained by 
applying Eq. d23l 



Wg° = h(u h -uJi) 



1 



_exp[P h hu h ] - 1 



exp[/3 ; frw;] - 1 

(57) 



C. Particle in an infinite square potential well 

For a QHE based on a particle confined in an Infinite Square 
(IS) potential, the eigen energies are E n (i) = JiU 2 , where 
7, = (Tih) 2 / (2mLf ); m and Li are the mass of the particle 
and the width of the square well at instant i, respectively. The 
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TABLE IV: Comparison of several properties of quantum Otto engines and quantum Carnot engines using the following quantum substances 
as working substances: Two-Level System (TLS), Harmonic Oscillator (HO), and a particle confined in an Infinite Square (IS) potential. The 
operation efficiency r\, the positive-working condition (PWC) and the amount of work W extracted in a thermodynamic cycle are listed below 
(see Sec. VI). 





Two-Level System (TLS) 


Harmonic Oscillator (HO) 


Infinite Square well (IS) 




Vo 


1 - — 




.-(ft)' 


Quantum Otto 
engine 


PWC 


T h > ^T t 


Th > %T t 


Th > T, 




Wo 




Wg° 


W Q S 




Vc 


1 - 

n 


1 - 


1 - 2L 

T h 


Quantum 
Carnot engine 


PWC 


T h > T t 


T h > Ti 


T h > T t 




W c 




Wg° 





entropy of the working substance can also be calculated by 
applying Eq. ( fT5T >: 



2 V 2 V PCii 

In obtaining Eq. (T58T > we have make an approximation 



^exp[ 7i n 2 ] 

71=1 



expj^n ]ofn. 



(58) 



(59) 



So the work done during a QCE cycle can be expressed as (fTTI i 



Wh s = (T h -T l )[S IS (B)-S IS (A)] 
while the net work done during a QOE cycle is 

-JS " 1 1 



Wo = o(7h-7i) 



{Phlh? (Pill) 2 



(60) 



(61) 



In order to better compare these results, these are all listed in 
Table IV. 

Before concluding this section, we would like to mention 
that, our current QHE model can only be implemented with 
quantum systems with all energy levels being discrete (all 
eigen states being bond states). Besides the harmonic oscil- 
lator and the infinite square well system, we can find other 
potentials that satisfy discrete spectral structure requirement 
to implement our QHE. However, we cannot deal with quan- 
tum systems with continuous spectral structure, e.g., a particle 
in a Coulumb potential or a finite square well. We will extend 
our current study to quantum systems with continuous spec- 
tral structure in future research. 



VII. MAXWELL'S DEMON AND QUANTUM OTTO 
ENGINE 



In the above discussions, we give clear definitions of mi- 
croscopic QCE and QOE cycles through clarifying the basic 
quantum thermodynamic processes (e.g., quantum isochoric 



process and quantum isothermal process). These results in- 
dicate that the properties, such as the operation efficiency, of 
macroscopic (classical) heat engines can be obtained from the 
microscopic (quantum) level as long as we clarify the basic 
thermodynamic processes microscopically. In the previous 
discussions, our QCE and QOE model show no contradiction 
with the thermodynamic laws. 

In one of our recent studies fisll . we proposed a Maxwell's 
demon assisted quantum thermodynamic cycle to study the 
function of a Maxwell's demon, and we also studied how it 
affects the second law of thermodynamics. It is interesting 
that when the restoration of the demon are properly included 
into the QHE cycle, the efficiency of the Maxwell's demon as- 
sisted QHE cycle has the same form as that for a QOE derived 
previously d27l ). Hence, the apparent violation of the second 
law due to Maxwell's demon is prevented. In this section, we 
also would like to study the intrinsic relation between these 
two kinds of quantum thermodynamic cycles. We would also 
like to add some details about the Maxwell's demon-assisted 
quantum thermodynamic cycle proposed in Ref. 11811 to better 
demonstrate our main idea. 



A. Maxwell's demon erasure not included in the 
thermodynamic cycle 

We first analyze a single-reservoir thermodynamic cycle 
with external control based on the effective temperature de- 
fined above. It can be proved that the property of this cycle is 
similar to that of the QHE cycle proposed in Ref. flal Our 
thermodynamic cycle consists of three steps: (1) quantum ro- 
tation, (2) decoherence, and (3) thermalization. 

After the thermalization process, the state of the two-level 
system (with the ground state |0), the excited state |1) and the 
level spacing A) can be described by a density matrix 



p(0) = Pi|l)(l| + P |0)<0|, 



(62) 



where the probability distributions Pi and Pq in the two-level 
system are determined by the temperature TJ of the heat reser- 
voir and the level spacing A. Then the two-level system, 
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driven by an external field (we do not treat it as a part of the 
system for the moment), undergoes a rotation, 



|1) -> |1) = cos0|l) +sin0|O), 
|0) -> |0) = -sin0|l> + cos0|O). 



(63) 



If the time interval of this rotation is much shorter than the 
relevant time of the thermalization, the state p(0) of the two- 
level system becomes 



p(l) = P{\1){1\ + P^|0)(0| + ODT, (64) 



where 



P[ = P 1 cos 2 9 + P sin 2 0, (65) 
P„ = Pi sin 2 6 + P cos 2 9, 

and ODT denotes the off-diagonal terms, which disappear 
rapidly (pure dephasing) due to the coupling between the two- 
level system and the reservoir. Actually two effects (dephas- 
ing and dissipation) occur when the two-level system is cou- 
pled to the reservoir [29]. The first process occurs much faster 
than the second one, and thus we can consider the two effects 
separately. After considering dephasing (but before consid- 
ering dissipation), the state p(l) of the two-level system be- 
comes 



p(2)=P 1 '|l)(l|+P^|0>(0| 



(66) 



This state is obviously not in equilibrium with respect to the 
reservoir at temperature Tj. But we can imagine there is such 
a reservoir at temperature Th, which can be expressed as 



T h (P ,Pi,t 



A 



In" 



Pi sin 2 6 + P cos 2 9 
Pi cos 2 6 + P sin 2 9 



(67) 



This effective temperature Th (Po , Pi , 0) possesses some ex- 
otic features. For example, when 9 = tt/2, P[ = P and 
Pq = Pi. This means that — — T) is a "negative tempera- 
ture" since there exists a population inversion P[ > Pq . Only 
when Pq > P[, Th is positive. Finally, the two-level system is 
put into contact with the heat bath for a sufficiently long time. 
After the thermalization process, the state p(2) returns to p(0), 
and a thermodynamic cycle is finished and the two-level sys- 
tem seems to extract work from a single heat bath. We can 
imagine it as a thermodynamic cycle between two reservoirs 
with the temperature Ti and a virtual temperature 2\. Actu- 
ally contradictions to the second law can appear due to this 
"negative temperature". 

It is not surprising that the above result (a contradiction to 
the second law) appears since we do not include the controller 
for the rotation operators shown in Eq. ( |63l . This result is 
very similar to those single-particle heat engines assisted by 
classical Maxwell's demon proposed by Szilard. 

Now, let us describe a new version of Szilard single-particle 
heat engine (see Fig. 8). The working substance of the QHE 
is a spin (two-level system) with ground state |0), and excited 
state |1). The level spacing is A. We first put the two-level 
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FIG. 8: (Color online) Schematic diagram of the two-level version of 
the Szilard single-particle heat engine. An invisible demon detects 
the state of the two level working substance and then controls it to do 
work: when the system is in the excited state, the demon makes the 
working substance flip to the ground state (the first step); when the 
system is in the ground state, the demon does nothing for the working 
substance. After these operations the working substance is brought 
in contact with a very large heat bath and then thermalized into its 
initial state (the second step). 



system into a heat bath at temperature T. After a sufficiently 
long thermalization process, the spin reaches a thermal equi- 
librium state, which can be described by p(0), similarly de- 
fined as above. Second, a demon performs a measurement. 
If the measurement result is confirmed that the system is in 
its upper state, then the spin is flipped and positive work is 
done by the spin with an amount A. If the measurement result 
confirms that the system is in its lower state, then no work is 
done. Then the system is put into contact with the heat bath 
and a new cycle starts. This is a two-step Maxwell's demon- 
assisted QHE. A similar discussion has been given in Ref. H. 
It is easy to see that the net effect of this QHE is to absorb 
heat from a single heat bath and convert it into work. On av- 
erage, the net work done per cycle is Pi A. This is a perpetual 
machine of the second kind. This apparent violation of the 
second law is seen because the erasure of the demon is not 
included into the QHE cycle. 



B. Maxwell's demon erasure included in the thermodynamic 

cycle 

In Ref. d, we have demonstrated that, when the erasure 
of the information stored in the Maxwell's demon is consid- 
ered into the QHE cycle, the apparent violation of the second 
law does not hold, i.e., there is no violation of the second law 
even in the existence of such a Maxwell's demon. 

To show the above observation, let us consider in more de- 
tail a thermodynamic cycle including the Maxwell's demon 
(see the Fig. 9). Let us explain each step of the QHE cycle 
proposed in Ref. lfl8ll . First, two qubits (two-level systems) S 
and D are decoupled and separately coupled to two heat baths 
with different temperatures T$ and To- After a period of time 
longer than both the dephasing time T2 and the relaxation time 




(b) (c) (d) 



FIG. 9: (Color online) Schematic illustration for the Maxwell's- 
demon-assisted quantum heat engine with four steps: (a): Both the 
QHE of the two-level system and a demon is hidden in the " piston" 
as another two-level system are initialized in thermal equilibrium 
states with different temperatures; (b) The demon makes a quantum 
non-demolition measurement about the state of the two-level sys- 
tem S with a CNOT operation: when seeing S in the excited state it 
records this information through a transition from its ground state to 
the excited state; when S is in its ground state, the demon remains 
in its original state; (c): According to the information recorded by 
the demon, the demon can enable the system to operate through an- 
other CNOT logical gate: If the demon was encoded in the excited 
state it will make the system flip; otherwise it enables the system to 
remain; (d) both the system and the demon are brought into contact 
with their own heat baths and then thermalized with different temper- 
atures from their own initial states, respectively. During this process, 
the information stored in the demon is totally erased and the entire 
system completes a quantum thermodynamic cycle. 



T±, they are thermalized to two equilibrium states pg(l) and 

Mi) 

Pf (1)=Pf\1)(M+Pf\0)(0\, (68) 

for F — S and D respectively. The joint thermalized state 
p(l) = ps(l) ® Pd(1) of the total system with S plus D can 
be written as 

P (l) = Pj; 1 |l,l)(l,l| + Pj;° |1,0}<1,0| (69) 
+P£ X |0,1)(0,1| + P££|0,0)(0,0|, 

where we have defined the direct product of the eigen state of 
the two qubits 

\q,q') = \q) S ®W)D, (<?,</ = 0, 1) (70) 

and the joint probabilities Pj'f, = Pf P^,. 

Second, two consecutive unitary operations: a CNOT op- 
eration flipping the demon states only when the working sub- 
stance system is in its excited state [28], and the demon con- 
trols the system to do work. Physically, the system experi- 
ences a conditional evolution (CEV). The changes of the states 
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in the two operations can be expressed as follows 

|0,0) CNOT |0,0) CEV |0,0) 

10,1) CNOT 10,1} CEV 10,1) 
' 1 ► 1 ' 1 > L' 1 , (71) 

1.0) CNOT 1 1,1) CEV |1,1) 

1.1) CNOT |1,0) CEV [1, 0) 

where |0) and |1) are defined in Eq. ((63). These two sub- 
processes can be realized with two quantum non-demolition 
Hamiltonians 

-Hcnot = Hs^D, Hqev = Hd^s, (72) 

where 

H A ^ B = 9 -{l + ai)®(j* (73) 

is a typical CEV Hamiltonian with A controlling B. It pro- 
duces a time evolution of the system B as a quantum rotation 
defined by Eq. d63l with 8 = gt, when the system A is ini- 
tially prepared in \1)a- Here, af, F = A,B, (j = x,y,z) 
are the Pauli matrices of the two level system A and B. Obvi- 
ously the CNOT operation on the demon is given by the time 
evolution produced by iicNOT with 9 = it/2 . iicNOT per- 
forms a quantum non-demolition measurement of the state of 
the working substance S. 

After this CNOT operation, the density matrix changes 
from p(l) to 

p(2) = CNOT{p(l)} (74) 

= P]; 1 |i,o><i,o| + pJ;° |1,1><1,1| 
+P°;J, |o,i)(o,i| + P°£|o,o)<o,o|. 

From Eq. d74b we can derive the reduced density matrices 
p s (2) and p^(2) of the working substance and the demon 

Ps (2) = Tr Z5 [p(2)]= Ps (l) (75) 

and 

p D (2) = Tr s [p(2)] (76) 

= (Ps'D+PsUVW + iP^D +P5,d)|1)(1| 

respectively. During this process, the system S remains in 
its initial state while the demon acquires information. Mean- 
while, the internal energy changes at the expense of using an 
amount of work Wd 

W D = Tx[Hp D (2)]-Tr[Hp D (l)} (77) 
= A^PA-P^-P"' 1 ). 

From Eqs. d69l and ( f74b we see that the total entropy of S and 
D does not change during the CNOT operation 

S(2) = -k B [Plf D hxP l s % + P^° D lnP 1 s: ° D (78) 

, p0,l i p 0,l , p 0,0 , p0,0i 

= -k B [Pi In P| + P° In P|] 

-k B [Ph^Ph + Pr>^PD] 
= S s (l) + S D (1) = S(1). 
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Hence, our model verifies the prediction in Refs. [33, 34] that 
a measurement do not necessarily lead to entropy increase. 
Also, we know from Eq. (l75l l that the entropy of S during the 
CNOT operation does not change 



S s (2) = S s (l). 
But the entropy of D changes as 

6S D = Sd(2)-S d (X) 



(79) 



(80) 



i,i 



= -k B [(P S 'D + PsId) 11 ^ 1 S.l> 



pO.O s 
r S,D) 



P 



0,1 



S,D 



HPs'°D 



We would like to point out that the mutual entropy 

S M (2) = S D (2) + S s (2) - S(2) 
= S D (2)-S D (1). 



(81) 



does not vanish (it vanishes before this CNOT operation). And 
this non-vanishing mutual entropy can be used to measure the 
information acquired by the demon about the system 1 30] . 

Next we consider the changes of entropy and energy after 
the quantum control (the CEV) process. The state of the de- 
mon and the working system after the CEV is 



p(3) = CEV{p(2)} 



Ps,d M) M + Pslh |M)(M| 
+Ps,d |o,i> <o,i| + |o,o) (o,o| 



(82) 
(83) 



The reduced density matrices of the working substance Ps(3) 
and the demon pp(3) are respectively 



p s (3) = Tv D [p(3)} 



(84) 



^|1>(1|+P^|0>(0| 
P°s'° D 10} <0| 



-Ps%\l)(l\ 



and 



p D (3)=Tr s [p(3)}=p D (2). 



(85) 



Similarly, the energy change (work performed by the working 
system) during this process is: 



W S = Tr[Hp s (2)] - Tr[H Ps (3)] 



(86) 



As ^ - P. 



S,D 



P^Ki|i)| 2 -P5 il<6|i>| 2 ). 

In particular, when we choose 8 — ir/2, the CEV is a CNOT, 
and the reduced density matrix of S can be written as 



Ps(3) = PA |i)(i|+p£ |0) <0|, 

and the entropy of the system is 

S s (3) = -k B [Ph toPh +PSlnP&] = S D (1) 



(87) 



So the change of the entropy of the working substance in the 
CEV process is 

SS S = S s (3) - S s (2) = S D (1) - S s (l), (88) 

where we have used Eq. ( l79b . But from Eqs. (F78T >. d82l and 
dHBT ). we know that both the total entropy of S plus D and the 
entropy of the demon D do not change, i.e. 



S(l) 



5(2) = S(3), 
Sd{2). 



(89) 



Finally, S plus D are decoupled and put into contact with 
their own baths, and a new cycle starts. In the thermalizaiton 
process, no work is done, but heat is exchanged between the 
heat baths and S and D. The thermalization process is an 
information-erasure process. This kind of zero-work erasure 
with a low temperature reservoir was first introduced in Ref. 
l3ltl . and studied afterwards in Ref. l32ll . In the thermaliza- 
tion process the heat absorbed by S is 



Q m = Tr[Hp s (l)}-Tr[Hp s (3)} 



(90) 



As(P 



1,0 

S,D 



l(i li>l 1<b li>l ). 



and the heat released by the demon D is 



Qo 



Tr[Hp D (3)]-Tr[Hp D (l)] 



(91) 



= An IP, 



,1,0 

S,D 



p0,l pi 



Now we include the erasure of the memory of Maxwell's 
demon into the thermodynamic cycle, and we will show that, 
under certain conditions, our composite QHE is equivalent to 
a simple QOE. 

For each cycle described above, we are now able to calcu- 
late the work W performed by the heat engine 



W 



Ws-W D 
As(P 



(92) 



,1,0 

S,D 



- P, 



1.0 



S,D 



(il 1 )! 



-p: 



S,D 



|(6|i>| 2 )-A D ( 



pi,u _ p 
r S,D r ' 



S,D 



It can be checked that in the thermodynamic cycle W = 
Qm — Qout- This is just the first law of thermodynamics. The 
positive-work condition can be derived from Eq. ( 1921 



T s > Tr 



(93) 



Notice that when we choose the CEV to be the special case 
9 = 7r/2, (i.e., a CNOT) the heat absorbed Q ln and positive 
work W done during a cycle can be simplified to 



Qu 



w 



(plfi 

\ S,L 

( P 1 ' 
\ S,L 



_ p0,l 
D r S,D I ' 



_ P 0,l 
D r S,D 



( P h ° 

ys.D 



p. 



S,D 



(94) 
(95) 



and 



W 



. ( pi.i _ pM" 

An VS.D r S,D, 

1- . ; (96) 
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If we further assume the temperature Td to be so low that 
exp(— P d Ad) <C 1, i.e., the demon is "erased" nearly to its 
ground state p D (1) « |0d) (0d| (The zero-entropy "standard 
state" was also discussed in Ref. ll3Til32j[33ll ). Then the effi- 
ciency of our QHE, Eq. d96i >, becomes 



r?=l- 



A, 



(97) 



This is the efficiency of a simple QOE cycle without 
Maxwell's demon, as shown in Eq. d27b . Otherwise the op- 
eration efficiency in Eq. d96l ) is less than the efficiency of a 
simple QOE cycle. This is because i) among all CEVs, the 
CNOT is the optimal operation to extract work, and ii) when 
Td is vanishingly small, the demon can be restored to a zero- 
entropy "standard state" to acquire information about the sys- 
tem in the most efficient way. 



Bose-Einsten condensation, of the working substance. We be- 
lieve these quantum effects will advance our understanding of 
the relation between quantum mechanics and thermodynamics 
and bring important insights into some fundamentle problems 
in quantum thermodynamics. 
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VIII. CONCLUSIONS AND REMARKS 

By defining the quantum version of basic thermodynamic 
processes, we study the basic properties of a QCE. To con- 
struct a QCE cycle, two preconditions about the quantum adi- 
abatic process are required: first, all energy gaps of the work- 
ing substance must change by the same ratio, such as in har- 
monic oscillators and infinite square well potentials. Second, 
the change of the ratio of the energy gaps must equal the ratio 
of the temperatures of the two heat baths. We find that the 
working efficiency for the QCE is the same as that of the clas- 
sical Carnot engine though the internal energy may change in 
the quantum isothermal process. We also study the properties 
of the QOE and compare these with the classical Otto engine, 
and we find that the efficiency and positive-work condition 
are the same when expressed in terms of temperatures (see 
Table III). Through comparing the thermodynamic cycles of 
the two QHEs, we clarify the relationship between them, and 
we demonstrate that the QCE (QOE) cycle can be modeled 
as an infinite number of small QOE (QCE) cycles. We also 
discuss some experimentally realizable physical systems that 
can be used to implement our QHE. Finally, through a gener- 
alized QOE, we demonstrate that there is no violation of the 
second law, even when there is a Maxwell's demon. This is a 
prediction of Landauer's principle ll3^[3~5ll36i[3~7i[38l l. 

Before concluding this paper, we want to emphasize three 
points. First, in our present study the working substance is al- 
ways assumed to be in a thermal equilibrium state, and quan- 
tum coherence is not considered here, our study is not related 
with the result that the efficiency of a QHE is less than its 
classical counterpart ?/q < ri c , as in Ref. lf28tl . Second, we 
only consider quasi-static process. Hence, the time intervals 
of these processes are infinitely long and the output power is 
vanishingly small. Recently, finite-time QHE cycles fl5l l39tl 
(nonzero output power) and friction-like behavior 1 16] of the 
QOE were studied, where the increase of power occurs at 
the expense of decreasing operation efficiency. Third, we 
will further extend our current study to QHE with quantum- 
many-body system as the working substance. In this extened 
case, we will consider the quantum statistical effects, e.g., the 



APPENDIX A: INVARIANCE UNDER ENERGY SHIFT 

The amount of positive work Wo done during a QOE cycle 
remains invariant under the uniform shift of all energy levels. 
If we shift all the energy levels E l n in Fig.l 

E l n = E l n - 5, (Al) 

the amount of work done during a QOE cycle becomes 

n 

where P l n are the occupation probabilities after the energy lev- 
els are shifted. It is easy to find that the occupation probabili- 
ties P l n remain invariant under such an energy shift 



Pi 



<E i n -m 



(A2) 



-i -l 



= p'. 



Thus, after the energy levels shift, the net work Wo done dur- 
ing a QOE cycle can be simplified to 



Wo = J2(P%-PM-E l n + 5) 

n 



(A3) 



which is just the net work Wo done by the working system 
during a QOE cycle before the level shift. Similarly we can 
prove that Wo remains invariant under the shift E% = E„ — S. 

Following the same way we can prove that other properties 
of the QHE, such as operation efficiency and positive work 
conditions, are invariant under uniform shift of all energy lev- 
els as well. Hence, we can simply assume the eigen energy of 
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the ground state to be zero Eq 



El 



0. This energy shift 



is convenient for our discussion about QHE later. This result 
can be easily generalized to any QCE cycle, because in Sec. 
V, we demonstrate that a QCE can be modelled as an infinite 
number of small QOE cycles. 



APPENDIX B: INTERNAL ENERGY OF THE WORKING 
SUBSTANCE IN QUANTUM ISOTHERMAL PROCESSES 



The internal energy of a system is 

U(t) = Tr[p(i)H(i)}. 



(Bl) 



In an isothermal expansion process, all energy levels change 
in the same ratio 

E n ^(E n , n = 0,l,2,.--; 

where £ is the ratio of energy levels and < ( < 1. When the 
energy levels change, the internal energy of the system can be 
rewritten as 

U (0 = E §77) exp[-/3 h CEn] , (B2) 

n V ^ / 



where 



Z{Q = Y,zM-fh<E n ] (B3) 



is the partition function. To test whether the internal energy of 
the system is invariant under the change of energy levels, we 
take the derivative 

= E^)( 1 -^^™)exp[-/3,C^] (B4) 



E„ 



For two-level systems, there is only one term in the sum 
over n, because we assume the eigen energy of the ground 
state to be zero. Thus Eq. ( IB4I ) can be simplified as 



dim 



E e 

z(0 



expb/^C^e 



1 



z(0 



(B5) 



where E e is the eigen energy of the excited state of the two- 
level system. The rhs of Eq. dB5l ) is obviously nonzero. 

For a harmonic oscillator with the eigenfrequency u>, we 
assume the eigen energy of the ground state to be zero (ne- 
glecting the ground state energy haj/2), then Eq. ( lB4t can be 
simplified to 



dU(Q 



tuv{i3 h c,hijj + 1) 



Tujf3 h ( 



(B6) 



exp^CM - 1 



The rhs of Eq. (IB6b is obviously nonzero. 

For a particle confined in an infinite square potential, E n = 
771 2 , where 7 has been defined in Eq. (1581 . Then Eq. (IB4I) can 
be simplified to 



dU(Q 



1^ W7\ exp[-P h C-yn } (B7) 



+0hC 



Z(0 

E ~z{Q e M-0h(jn 2 } 



2 



where 



z(o-E ex p[-^C7n 2 ] 



(B8) 



is the partition function. We make an approximation 

f°° 1 I 

Z(Q « I eM-PhCin 2 ]dt = 2^/7^. (B9) 



and then Eq. dB7b can be simplified to 



dU(Q 
d( 



1 



(BIO) 



The rhs of Eq. ( IB 101 ) is nonzero, either. Hence, we conclude 
that in an isothermal process, the derivation of the internal 
energy over the energy levels in Eq. (IB4I ) is always nonzero. 
Accordingly, the internal energy of the system varies with the 
change of the energy levels. 



APPENDIX C: THERMODYNAMIC REVERSIBILITY OF 
AN INFINITE NUMBER OF INFINITESIMAL QUANTUM 
ISOTHERMAL PROCESSES 

We consider a two-level system interacting with a heat bath. 
The temperature of the heat bath is well controlled so that the 
two-level system is always in thermal equilibrium with the 
heat bath. Now let us calculate the entropy (Von Neumann en- 
tropy) increase in the two-level system and the entropy (ther- 
modynamic entropy) decrease in the heat bath. We assume the 
two-level system initially in thermal equilibrium with a heat 
bath at the temperature (see Eq. ( |40b and Fig. 6), and fi- 
nally the temperature of the heat bath is controlled to increase 
to Th (see Fig. 6). The entropy of the initial state of the system 
at A' (see Fig. 6) is 

Stls(A') (CI) 
= -k B (A'_ expI-Afc/^] m(A'_ expf-A^) + A'_ In A'_) 

= ^AV-fesMA'+exptA^]), 



where /3' h = l/(k B T^) and 



A' 



1 



A± 



1 + exp[±A h p' h ] 
1 

1 + exp[±Aft/? h ] ' 



(C2) 
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Similarly, we obtain the entropy of the final state of the two- 
level system at B' (see Fig. 6) is 

Stls(B') = - k B ln(A+ cxpIA/AD- (C3) 

J-h 

So the entropy increase of the two-level system in the in- 
finite number of infinitesimal quantum isothermal processes 
{A 1 — > B' in Fig. 6) is 



We apply the transformation 



T 



k B In t ' 



(C6) 



dT = *L 
k B 



ln 2 t 



dt 

T' 



Under this transformation dSeath dC5l > can be further given as 



^Stls — Stls(B') — Stls(A') 



(C4) 



-f A+~k B ln(A + exp[A^]) 
J-h 

^A' + -fc B ln(A' + exp[A^]) 

L h 



Next we calculate the entropy decrease of the heat bath in the 
same process. The thermodynamic entropy decrease of the 
heat bath due to its coupling to the two-level system can be 
expressed as 



dS- 



Bath 



h &hdpe 



1 1 i T 

± h 

k B / 

./Tl 



A^ 2 lk , 2 fA h ^ 2 



(C5) 

dT 
~T' 



dSB 



ath 



(C7) 



lni 

T+t 



exp[A h p h ] 
exp[A fc /3y 



hi 





exp[A h /3 h ]~ 




cxp[A h /3' h ]_ 



= -^A + - k B In (A+ cxp[A h /? h ]) 

J-h 



yA' + -k B lnAVcxp[A^] 



From Eqs. (1C41 > and dC7t we can see that dSsath = ^Stls* 
i.e., the entropy decrease in the bath is equal the entropy in- 
crease in the two-level system. Thus we proved the total en- 
tropy conserves in these infinite number of infinitesimal quan- 
tum isothermal processes, and these processes are thermody- 
namically reversible. 
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